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We formulate scattering in one dimension due to the cou- 
pled Schrodinger equation in terms of the S matrix, the uni- 
tarity of which leads to constraints on the scattering ampli- 
tudes. Levinson's theorem is seen to have the form 77(0) = 
7f(ni, + |n — ^N), where 7/(0) is the phase of the S matrix 
at zero energy, rib the number of bound states with nonzero 
binding energy, n the number of half-bound states, and N the 
number of coupled equations. In view of the effects due to the 
half-bound states, the threshold behaviour of the scattering 
amplitudes is investigated in general, and is also illustrated 
by means of particular potential models. 



I. INTRODUCTION 

The quantum mechanics of one-dimensional scattering 
describes many actual physical phenomena to a good ap- 
proximation. (For example, see Ref. for a review of 
tunneling times.) One-dimensional models are further- 
more often employed to gain deeper insight into the ap- 
proximations used in order to make the more complex 
three-dimensional systems tractable. It is therefore not 
surprising that there have been many articles, also in this 
journal, dealing with various aspects of such scattering. 
In particular a number of papers have appeared in recent 
years on the threshold behaviour of one-dimensional scat- 
tering and Levinson's theorem J^-Q]. These studies have 
been limited to systems without coupling. 

In this paper we wish to investigate scattering de- 
scribed by a system of coupled differential equations 
with a particular interest in developing a formulation 
for Levinson's theorem and in gaining insight into the 
threshold properties of the scattering amplitudes. This 
work can be seen as a special case of multichannel scat- 
tering for which the threshold energies are equal. In sub- 
sequent work, we intend to generalize to the case of dif- 
fering threshold energies. Although in previous work on 
one-dimensional scattering one has at times employed a 
"partial wave" analysis |t],|| or a parity-eigenstate rep- 
resentation [EMI, we have chosen to use the traditional, 



more "physical" , approach involving states with incident 
waves coming from a single direction. 

In Sec. we express the scattering properties in terms 
of the S matrix, the unitarity of which leads to speci- 
fiable constraints on the scattering amplitudes. For the 
proof of the generalized Levinson's theorem we make use 
of the complete set of orthonormal states of the Hamil- 
tonian; this is an alternative to the approach involving 
the analyticity of the scattering amplitudes. The proof 
of Levinson's theorem depends on the threshold proper- 
ties of the scattering amplitudes. These properties are of 
interest in their own right, especially in connection with 
scattering time delay and advance Q, and therefore we 
discuss the zero-energy behaviour of the amplitudes at 
some length. 

The factorization of the S matrix is generalized to the 
coupled system in Sec. III. We also indicate that there is 
a class of finitely periodic matrix potentials for which the 
scattering amplitudes can be found in a way analogous 
to the case with no coupling. 

In Sec. IV we discuss a number of specific potential 



models to elucidate and amplify general results. We con- 
clude with a brief discussion of our results in Sec. |V|. 



II. S-MATRIX FORMULATION 

The one-dimensional scattering problem has been 
studied in terms of the S matrix by a number of au- 
thors. (See, for example, Refs. [7p|-p^|.) We extend the 
formalism to include a matrix potential function. The 
Schrodinger equation for a stationary state of such a sys- 
tem is 



dx 2 



(2.1) 



where V(x) is a real, symmetric Nx N matrix, k 2 the en- 
ergy of the system, and ^(fc, x) the wave function, which 
is an ./V-dimensional column vector. For large values of 
I a; I the potential matrix V(x) approaches zero sufficiently 
fast, so that in the asymptotic region ^(k,x) represents 
a free-particle wave function. To ensure this we will take 



1 



for \x\ > R, R being the range of the poten- 



we assume that |V^'(a:)| is integrable 



V(x] = U tor \x 
tia(j. Furthermore 
for i,j = 1,...,N. 

The physical scattering solutions of Eq. (2.1) at a given 
energy k 2 can be written as the columns of TV x TV ma- 
trices ip(k,x) and ip(k,x) which are uniquely determined 
by the boundary conditions, 



and 



ip{k,x) 



r(k)e 



tkx 



X — > oo 



X — > — oo 



+ p(k)e lkx , x — > oo. 



(2.3) 



We will refer to -0 and V" as t ne solution matrices. Note 
that the columns of ip contain the wave functions with 
an incident wave from the left, whereas the columns of 
if) includes those with an incident wave from the right. 
The TV x TV matrices p, p, r, r are generalizations of the 
usual reflection and transmission amplitudes [ ^3|Jl4| . The 
set of TV-dimensional column vectors of matrices ip and ip 



represent solutions of Eq. (2.1). The linear independency 
of these solutions can be shown by considering the 27V x 
2TV matrix, 



w(ip,ip) 



ip if) 
ip' ft 



(2.4) 



in which the prime indicates the derivative with respect 
to x. In Appendix A we show that the det W(ip, ip) is 
a constant, which is nonzero if and only if the solutions 
are linearly independent. Its value, determined from the 
asymptotic forms of ip and ip, is 



detW(ip,ip) = det(r(fc))(-2ifc) 



N 



(2.5) 



Thus when both k and det r(fc) are nonzero, the columns 
of ip and ip give 2TV linearly independent solutions. 

In order to define the S matrix, we consider the general 
solution matrix in the asymptotic region, 



Ae lkx + B'e- lkx , x 
A'e lkx + Be~ lkx , x 



oo, 



(2.6) 



where A,A',B,B' are N x N matrices. Since the un- 
primed matrices are associated with incoming waves and 
the primed matrices with outgoing waves, the S matrix 
can be defined as the matrix that transforms the coeffi- 
cients of the incoming waves into those of the outgoing 
waves so that 



1 The boldface and 1 refer to the zero and identity N x TV 
matrices, respectively. 



A' 
B' 



= S 



Sba Sbb 



A 
B 



(2.7) 



Clearly, S is a 27V x 2N matrix. We can write it in terms 
of the transmission and reflection amplitudes by making 
use of the special cases for which (A, B) = (1,0) and 
(A,B) = (0,1). The result is 



S = 



( S a a 


S a b \ 




T p \ 


\ Sba 


Sbb J 


"( 


v P 7 J 



(2.8) 



The S matrix contains AN 2 complex elements or 8iV 2 real 
parameters. As we will see below there are a number of 
relations between the transmission and reflection ampli- 
tudes, which will reduce the number of independent real 
parameters. 



A. Relation between reflection and transmission 
amplitudes 

Constraints on the transmission and reflec tion ampli- 
tudes follow from the Schrodinger equation (|]l]). Con- 
sider two solution matrices ip±(k,x) and ip2(k,x), then 

ip%=(k 2 -V)ip 2 and 4" =iPl(k 2 -V), (2.9) 
so that 

4"iP 1 = 4(k 2 -V)ip 1 . (2.10) 



The Schrodinger equation, Eq. ( |2.l| ), for tpi premulti- 
plied by ip\ yields 



iP\iP'{ = iPl{k 2 -V)^. 
Subtracting Eqs. ( |2.10| ) and fl2.1l| ), we obtain 

4'^ - 4ip[ 



(2.11) 



- 441 



d 
dx 



which leads to 



4'i J i ~ 4^1 = constant matrix. 



0, (2.12) 



(2.13) 



If we now inser t the asymptotic forms of ip or ip for ip\ 
or ip2 into Eq. ( [2.13 ) and equate the expression at — oo 
to that at +oo, we obtain the following relations. 



t't + p^ p = t^t + pi p= 1, 



(yr + p - 



p + jCMT= 



(2.14) 
(2.15) 



Eqs. (2.14) and ( ^.15 ) are equivalent to the statement 
that S T S = I, where I is the 27V x 27V identity matrix. 

Further relations are found by using the time-reversal 
symmetry of the system. Since V(x) is real, the com- 
plex conjugate solution matrices ip* and ip* are also so- 
lutions of the Schrodinger equation. By complex conju- 
gating Eq. ( |2.6| ), we see that the roles of the incoming 
and outgoing asymptotic waves are reversed, and indeed 
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Eq. <^M is valid when A — > B'*, B — > A'*, A' 
and B' —> ^4* . Thus Eq. (2/7) may be written as 



.4* 



5 



A'* 



B* 



(2.16) 



Multiplying on the left by S> and using S*S — I, we find 
that 



B'* 
A'* 



which leads to 

A' 
B' 



1 





s 1 



B* 
A* 



1 





where S T is the transpose of S. Thus 



S = 



1 





s 1 



from which it follows that 



and 



1 





(2.17) 



(2.18) 



(2.19) 



(2.20) 



Inserting these expressions into Eqs. ( ^.14| ) and ( [2.15 ) 
we obtain 



rp' 



pr 1 = pr 



W= i, 

+ Tp= 0. 



(2.21) 
(2.22) 



These equations yield the other half of the unitarity con- 
dition of the S matrix, so that 



S*S = SS^ = I. 



(2.23) 



For a parity-invariant potential function, i.e., V(—x) = 
V(x), there are further constraints on the transmission 
and reflection amplitudes. In that case the amplitudes 
are symmetric matrices and the two types of amplitudes 
are the same, i.e., 



~ T ~T 
P=P=P =P , 



and 



~ ~T T 
T = T = T = T . 



(2.24) 



There is also a useful relation between the scatter- 
ing amplitudes at k and at — k, which is easily obtained 
by generalizing the result for uncoupled potentials 0. 

Since ip*(—k,x) and tp*(—k,x) are solution matrices of 
the Schrodingcr equation with the same boundary con- 
ditions as il)(k,x) and tp(k,x) respectively, we find that 



r*(-fc) =r(k), p*(-k)=p(k), 
r*(-k)=r(k), p*(-k)=p(k)- 



(2.25) 



From these relations it follows immediately that the re- 
flection and transmission amplitudes at threshold (i.e., 
k = 0) are real. 



B. Levinson's Theorem 

Levinson's theorem in its most common formulation 
for a spherically symmetric potential gives a relationship 
of the scattering phase shifts at zero and infinite energy. 
The theorem has also been studied for one-dimensional 
systems without coupling We generalize the the- 

orem to the matrix-potential case. Levinson's theorem 
is a consequence of the orthogonality and completeness 
relation of the eigenstates of the total Hamiltonian. 



The scattering state s of th e Sc hrodinger equation (2.1) 
defined by Eqs. ( |2.2| ) and (2.3) along with the bound 
states can be used to form a complete orthonormal set 
of eigenstates. Suppose that there are rib bound states 
whose orthonormal wave functions are denoted by the 
column vectors ipEj{x) with Ej(< 0) referring to the 
bound state energy. In case of degenerate bound states, 
we label the state vectors with subscript Ej , but allow the 
possibility of different subscripts j for the same energy in 
order to include all independent bound-state vectors. For 
example, one could have Ej = Ei where i ^ j when Ei 
is a degenerate energy eigenvalue. The orthonormality 
relations are 



2tt 



dx ipj(k, x)i/Ji(k', x) 
1 f°° ~ ~ 

— j dxi>}{k,x)i) i (k',x)=6 ij 5(k' -k), (2.26) 



/2ttJ_ 



dx tpUk,x)tp Ei (x) 



2tt 



dx ip](k,x)ipEi(x) = 0, 



(2.27) 



1 f°° ~ 

— I dx ip](k, x)ipi(k', x) — and 

27T 7 ^ 



dx ^ E .(x)ipEi(x) = Si 



(2.28) 



where ipi and ipi are the i th columns of the ip and i/j 
matrices respectively. Thus the completeness relation is 

5>e 3 (z)<V) + E^- / dk^{k,x)4{k,x') 

j=l i=l 77 Ja 

N 1 r°° 

+ V — / dktpi(k,x)tpUk,x') = lS(x-x'). 

— ' 27T In 

i=l J{ > 

(2.29) 

The completeness relation may be written in a more com- 
pact form in terms of the matrices themselves rather than 
the column vectors, i.e., 

V ip Ej . (x') + — / dk [V>(fc, x)^ (k, x') 

1=1 2n J ° 

+ ^j(k,x)ft(k,x')} = lS(x-x'). (2.30) 



3 



For the free-particle case when V(x) = 0, there are no 
bound states and the completeness relation is 



1 

2^ 



dk[ip°{k, x)i/j oi {k, x 1 ) + tp°(k, x)ip°Hk, x 1 )} 
= lS(x-x'). 



(2.31) 



and 



21 dk = Tr 



= Tr 



diA, 



dk 



= Tr 



of 



= Tr 



5 dk 



— t t H — -p T H t t H — -p 1 

9fc 9fc <9fc <9fc 



(2.38) 



We now subtract Eq. (|2.3l|) from Eq. fl2.30|) , set cc' = 
x, and integrate over x from —72. to 72. The resulting 
equation may be written as 



dk / dx [ip(k, x)ijj' [k, x) + ip(k, x)ip'{k, x) 

-i>°(k,x)^(k,x) -^°(k,x)^(k,x)} 

/n "b 
dasV^CajjVi.C*)- ( 2 - 32 ) 



The trace of Eq. ( |2.32| ) in the limit as 72 approaches in- 
finity gives 



lim / dk I i 

R ^°°io J-TZ 



dx Tr[^(fc, x)ip^(k, x) + ip(k, x)ip^{k, x) 



-ip (k,x)ip° J '(k,x)-^ (k,x)ip^(k,x)] =-2im h (2.33) 

To p erfor m the integration over x in the right side of 
Eq. ( 2.33 ), we use the identity 



dip &ip^ d 2 tp | 
dk dx dxdk 



(2.34) 



which may be obtained by taking the derivative with re- 
spect to k of the Schrodinger equation fl2.ip . Since in the 
limit 72 exceeds the range of the potential, we can insert 
the asymptotic forms of the wave functions in Eq. ( 2.33| ) 
to obtain 



Thus Eq. (2.35) may be written as 



r;(0) — 77(00) = nut — lim X(TZ), 

1Z — >oo 

where 

1 f°° dk 

-i(pHk) + p{k))e- 2lkn 
dk 

: Tv[p(k)+p{k)\e 



(2.39) 



2ikTZ 



2ikTZ 



k 



(2.40) 
(2.41) 



In the next section we show that p(k ) ~ Q (l/k) for large 
|fc|, so that the integration in Eq. ( 2.41 ) converges for 
large |fc|. We now take the limit as 72 approaches 00, 
using the relation 



^2ikTZ 



lim — — 

TZ^oo k 



= inS(k), 



(2.42) 



where S(k) is the Dirac delta function. In that limit X(TZ) 
goes to — (7r/4)Tr[p(0) + p(0)], so that the statement of 
Levinson's theorem now is 



r ? (0)= 7 rn 6 + -Tr[p(0)+p(0)], 



(2.43) 



where we have set 17(00) equal to zero. In the next section 
we also show that 



Tr[p(0)+p(0)] = -2(iV-n) ) 



(2.44) 



f°° dk m 

hm / — Tr 

iz^oo J 2k 

2ik1Z 



i(p + p)e zlKK +z(p T +p T )e 



' '[jJk T 
t , rt\ 



dr + dp + 9t_ + dp^ 
dV dk dk p 

= -2nn b . (2.35) 



-2ikTZ 



where n is the number of "half-bound states" |10] , |1 1| , |18| ] 
Thus in its final form Levinson's theorem states, 



1 



1 



77(0) = ir(n b + -n- -N) 



(2.45) 



Following Newton and Jost m& we define the phase as 



T)(k) = ^-lndctS(fe), 
2i 



(2.36) 



where we require rj(k) to be continuous for k G (0, 00). 
Since the 5* matrix is unitary, we may write S = U' SdU 
where U is a real orthogonal matrix and Sd is the diag- 
onal matrix Sn = diag(e j51 , . . . , e lS2N ) where the <5j's are 
real phases |l7|| . Let us write therefore S = e lA where 
A = WA D U for A D = diag(<5i, . . . , 6 2N ). Then 



2ir)(k) = lndetS(fc) = In det S D (k) 



2N 



j'=i 



(2.37) 



This expression of the theorem is consistent with that for 
the uncoupled case given in Ref. flT||. 



C. Threshold behaviour of p and r 

The threshold behaviour of reflection and transmission 
amplitudes and coefficients has been discussed recently in 
several articles . In order to study the behaviour 

of the p and r matrices at k = 0, we introduce a dif- 
ferent set of solutions of the Schrodinger equation, since 



according to Eq. (2.5) the columns of ijj and ip fail to be 
linearly independent at k — 0. Let </>(fc, x) and x(^> x ) be 
solution matrices of Eq. (2.1) which satisfy the boundary 
conditions, 
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<t>(k-R) = X '(k,-R) =1 
<f>'(k,-R) = X (k,-R) =0, 



(2.46) 



where R is the range of the potential as defined at the 
beginning of Sec. II. By evaluating det W(tf>, x) at % — 
—R and using the results of Appendix A, we readily show 
that for all x and k the de t W (<t>, x) = 1, where the matrix 
W is defined as in Eq. (^4). Thus unlike the column 
vectors of ip and the column vectors of <j) and x are 
linearly independent at zero energy. 

In order to obtain the scattering amplitudes, we ex- 
pand ij) in terms of <f> and x so that 



ip(k, x) = <j>(k, x)B{k) + X (k, x)C(k), 



(2.47) 



where B(k) and C{k) are matrices of expansion 
coeffic ient s. Evalu ating ip and tp' at ±i? using 
Eqs. ( |2.2|) and ( 2.47 ), we obtain four equations involv- 
ing B (k) , C '(k) , p(k) and r(k). By eliminating three of 
these we find that p can expressed in terms of <p ancl X 
and their derivatives evaluated at R. Thus 

p(k) = {k 2 x (k, R) + ik[ x \k, R) + R)] - cf>'(k, R)}- 1 
{k 2 X (k, R) + ik[ X '(k, R) - <j>(k, R)] + 4>'(k, R)}e~ 2lkR . 

(2.48) 

Similarly by expanding ij> in terms of <f> and x we obtain 
p(k) = {k 2 X (k, R) - ik[ X '(k, R) - (f>(k, R)} + <P'(k, R)} 

-l-2ikR 



{k 2 X (k, R) + ik[ X '(k, R) + 4>(k, R)} - 4>'{k, R)} 



(2.49) 



r(Jfc) = 2ik{k 2 x (k, R) + ik[x'(k, R) + 4>{k, R)] 

~4 > '{k,R)y 1 e' 2lkR , (2.50) 

and, since r(fc) = ^(k), 

r(k) = 2ik{k 2 x T (k, R) + ik[x' T (k, R) + (f> T (k, R)} 

-<t/ T (k,R)}- 1 e-* kR . (2.51) 



Thus all four scattering amplitudes, and consequently the 
S matrix, are determined by <f> and x an d their deriva- 
tives evaluated at R. In Appendix B it is shown that 
(j)(k,x) and x{k,x) are entire functions of complex k for 
all x e [-R, R] so that the analytic properties of the scat- 
tering amplitudes can be determined using these wave 
functions. By inserting the expressions for th e lar ge rea l 
k behaviour of the wave functions, Eqs. (Bl£) and ( B17 ), 
into the expressions for the scattering amplitudes, we ob- 
tain 



r(k) 



0(1/ k) andp(fc)-0 



0(1/ k) for k -> oo, 
(2.52) 



and similar expressions for r and p. 

If det[<^'(0, R)] 7^ 0, the reflection and transmission 
amplitudes at zero energy are 



p(0) = p(0) 



and t(0) = r(Q) 



(2.53) 



The case for which det[0'(O, R)] — needs special at- 
tention. In order to understand the significance of this 
condition, we look at the b oun d states of the system. In 
the Schrodinger equation ( fO] ) for bound states, we de- 
note the bound-state energy as a 2 — —k 2 with a > 0. 
The bound-state wave functions can be expressed as col- 
umn vectors of a matrix ipi,(a,x) with the asymptotic 
boundary conditions, 



ipb(a, x) 



x < -R 



(2.54) 



where Q and T are matrices of constants. The number of 
independent bound states at a given energy will depend 
on the rank of the matrix tpb, and consequently cannot 
exceed N. Proceeding as we did for the scattering states, 
we expand the bound-state wave functions in terms of 
the functions <pb(a,x) and X&(a, x) which are solution 
matrices of the Schrodinger equation with energy —a 2 
and satisfy the boundary conditions 



4> b (a,-R) = X ' b (a,-R) =1 
Xb(a,-R) =ct>' b (a,-R) = 0. 



(2.55) 



Thus 



ipb(a, x) = (j)b(a, x)P(a) + Xb(a, x)j(a), (2.56) 

where f3 and 7 are matrices of expansion coefficients. At 
R and — R we mat ch the asymptotic form of the wave 
function, E q. ([2.54 ), and its derivative to the expanded 
form, Eq. ( 2.56j ), and its derivative. Eliminating T, Q 
and 7 from the four equations so obtained, we are left 
with the equation 



{a 2 Xb(a, R) + a[ X ' b (a, R) + <f> b (a, R)] 



c/>> b (a,R)}p(a) 







(2.57) 



Since one of the four matching equations is Q = (3(a)e aR , 
there will be bound states only if the matrix (3(a) con- 
tains nonzero entries. Such a nontrivial matrix exists 
only when 



det{a 2 Xb(a, R) + a[x' b (a, R) + (j> b (a, R)] 
+ cf>' b (a,R)}=0. 



(2.58) 



This is the bound-state eigenvalue equation for energy 
—a 2 . In contrast to the case with no coupling, the bound- 
state eigenenergies can be degenerate. 

Let us consider the eigenstates when a = 0. These 
will occur only if det[4>' b (0 , R)} = 0. In general the solu- 
tions represented by the columns vectors of ^(0, x) are 
bounded but not square integrable; hence they are re- 
ferred to as half-bound states |l^] . The restriction on the 
potential function that it vanishes for |a;| > R, precludes 
the possibility of having normalizable state functions at 
zero energy. For this to be the case a linear combination 



5 



of the columns of ipb would yield ^(0, x) = for |ir| > R. 
Such a bou ndar y condition would lead to the trivial solu- 
tion of Eq. (2.1). Normalizable zero-energy bound states 
can exist for potentials which are less restrictive than 
those of this paper [Q. 

Since 0(0, x) and </>&(0, x) are solutions of the same 
system of differential equations and both have the same 
boundary conditions, <fib(0, x) = <fi(0,x). Thus the condi- 
tion that 



det#'(0,i?)] = 



(2.59) 



is equivalent to the condition for the existence of half- 
bound states. 

Consider the matrix eigenvalue equation, 



<f>'(0,R)P = A/3, 



(2.60) 



where j3 is a column vector and A is its eigenvalue. There 
will be a nontrivial solution only if 



det[0'(O,i?) - Al] = 0. 



(2.61) 



Suppose the eigenvalues obtained are Ai,---,Ai\r- At 
least one of these must be zero if there is a half- 
bound state. Actually there may be n(< N) zero 
eigenvalues. We can order these in the following way: 
0, . . . , 0, A„+i, . . . , Xn- These n zero eigenvalues will 
have n linearly independent eigenvectors associated with 
them, which represent n distinct half-bound states. 
We now return to the discussion of reflection and trans- 



mission amplitudes. The inverse of r of Eq. ( 2.5C ) can 
be written as 

2ikr- 1 (k) = {k 2 X (k, R) + ik[ X '(k, R) + (j>(k, R)] 

~<P'(k,R)}e 2lkR . (2.62) 

If this equation is combined with Eq. ( 2.49Q , we obtain 

2ikp{k)T~ 1 {k) = {k 2 x (k, R) - ik[x'(k, R) ~ &(k, R)} 

+4>'(k,R)}. (2.63) 



The factor in the curly brackets of Eq . ( 2.62| ) is precisely 



that in the determinant of Eq. (2.58) (with k=ia), i.e 



it is the factor which determines the bound states of the 
system. In Appendix B we show that matrices <fi and x 
are elementwise entire functions of k. From Eq. ( 2.62| ) 
and the fact that r T = r we see that 2ikr~ 1 (k) and 
(2ik) N detr~ 1 (fc) are als o ent ire functions of k. Accord- 
ing to Eqs. ( 2.5S ) and ( 2.62 ) detr _1 (fc) has a zero at 
k = ia when —a 1 is the energy of a bound state. Thus 
det r(k) has poles (possibly multiple) on the positive 
imaginary axis of the complex k plane. In the absence of 
a half-bound state, detr _1 (fc) has an N th order pole at 
k = and detr has an ./V th order zero at k = 0. Since 
there can be no more than N half-bound states det r(k) 
and r(k) are analytic in a neighborhood of k = 0. 

In the following we con sider real k . Ta king the limits 
as k goes to zero of Eqs. ( 2.62) ) and (2.63) we obtain 



lira 2ikr- 1 (k) = -<j>'(Q,R) and 



fc— o 



lim 2ikp(k)r~ 1 (k) = (f>'(0,R). 

k — *0 



(2.64) 



Similarly using Eq. ( 2.51 ) and the transpose of Eq. ( 2.48 ) 
we get 



lim 2ifcr~ 1 (fc) 



(0, R) and 



lim 2ikp(k)T~ L (k) = ft 1 (0,R). 



(2.65) 



We introduce unitary matrices U(k) and V(k) which di- 
agonalize r [^o], page 192], so that 



TD (k) = U\k)T(k)V{k) and 
2ikr D 1 = V\k)2ik T - 1 (k)U{k). 



(2.66) 



In the limit as k approaches zero U(fi) and V(0) also di- 
agonalize </>' T (0, R), i.e., 0g"(O, R) = ^(0)^(0, R)U(0), 
so that 



lim 2ikTp\k) = -(j)' D T (0,R). 

k — >0 

We define the matrix 

r(Jfc) = V f (k)p(k)V(k), 



so that 



lim 2ikr(k)T D - 1 (k) = 4>' D (0,R). 



fe^O 



Combining Eqs. ( ^.67| ) and (2.69) gives 



r(0W D T (0,R) = -<f>' D T (0,R). 



(2.67) 



(2.68) 



(2.69) 



(2.70) 



The matrices <p' T (0,R) and <f>' D T (0,R) have the same 
rank (2(], page 55]. Thus cj)' D T (0,R) will have the same 
number of nonzero diagonal elements as there are nonzero 
eigenvalues of <j)' T {0, R) or 0'(O, R). Writing f n T (0, R) = 
diag(0, • • • , 0, s n +i, • • • , sn) and using Eq. ( 2.70 ) we see 
that the matrix r(0) must have the form 



r(0) 



R u 
R 21 -1 



(2.71) 



where the matrices Ru, R21, and 1 have dimensions 
n x n, (N — n) x n, n x (N — n) and (N — n) x (jV — n) 
respectively. 

To study the behaviour of ro(fc) near k — 0, we 
consider the Hermitian positive semi-definite matrix 
T(k) = T^(k)r(k), whose real nonnegative eigenval- 
ues we denote by t\ (fe), ... , t 2 N {k). The singular values 
of r(k) are defined as the nonnegative square root of 
these, i.e., . . . ,tisi{k), and they form the diago- 

nal elements of To(k) in Eq. ( 2.66| ), so that ro(fc) = 
diag(fi(fc), . . . ,tjs[(k)) H^, page 192] . Since r(fc) is ana- 
lytic in the neighborhood of k = 0, so is T(k). For small 
real k > 0, we invoke a theorem of Rellich |2l], page 
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31] which states that the eigenvalues of T(k) are c onver - 
gent power series of k. Furthermore, using Eq. ( 2.25 ), 
we note that T(k) — T T (—k) and the eigenvalues of 
T(k) and T(—k) are the same. Thus taking the eigen- 
values in the same order, we find that tf(— k) — tf(k). 
Hence the t 2 (k)'s are power series of k 2 and the tj(fc)'s 
are power series of k. Thus we may write Tn(k) — 

diag(£i + tnk + •••,..., t N + t N1 k H ) and T^ j 1 (k) — 

diag((ii + tnk H . . . , (t N + t N ik H where 

the ti's and ty's are constants. In order that 

lim 2ifcr £) 1 (fc) = lim diag(2ifc(ii + tnk H )~\ . . . , 

k^O k—*0 

2ik(t N +t N1 k H (2.72) 

= diag(0, . . . , 0, -s„+i, . . . , -Sjv), (2.73) 

the quantities t\,...,t n ^ and t n +\ = ■■■ = tN — 
0. The matrix T£>(0) therefore has the form 175 (0) = 
diag(ii, . . . ,t n , 0, . . . , 0) where the first n diagonal ele- 
ments are nonzero. 

Recall that t t = r and therefore ( r 7 ") - 1 =t~ 1 , and 
that p = 'p r . The second part of Eq. (2.64) becomes 



lim 2ikr' 1 (k)p{k) = /T (O, R). 

k^Q 

The same U and V can be used to obtain 
lim 2ikT D - 1 (k)r(k) = ^ T (0,i?), 



fc^O 



where 



r(k) = U\k)p(k)U{k). 



(2.74) 



(2.75) 



(2.76) 



As before the structure of r(0) may be determined, and 
it is 



r (0) = 



-1 



(2.77) 



In order to simplify Eq. ( 2.43| ) we write 

Tt[ P {0) + p(0)} = Tr[r(0) + r(0)] 

= -2(A^-n)+Tr[i?' 1 \ x ™ 



i?™^"], (2.78) 



where the superscripts refer to the dimensions of the ma- 
trices. All that remains is to evaluate the last term of 
the right side of this equation. At zero energy the scat- 
tering amplitudes are real. Using this fact along with the 
relations that r = r T and p = p T in Eq. ( [Of ), we obtain 



-(0)p(0)+p(0)T(0) = 



(2.79) 



By applying the diagonalization transformation of r, we 
find that 

r D (0)r(0)+f(0)r D (0) =0. (2.80) 
In matrix form this equation may be written 



Ai W i? n 
U21 -1 



i?H i?i2 
1 







(2.81) 



where A\ = diag(^i, • • • , t n ) and the bloc k ma trices have 
the appropriate dimensions. From Eq. (2.81) it follows 
that 



Aii?n +R n Ai = r 



(2.82) 



Cons eque ntly, the second term in the right side of 
Eq. (2.78) is zero, and we have the simple relation 



Tr[p(0) + p(0)] =-2{N-n) 



(2.83) 



In light of the discussion of the N = 1 normal and 
anomalous threshold behaviour HJtI], our result seems 



surprising since the right side of Eq. ( 2.83 ) is an inte- 
ger, whereas the anomalous threshold effect for parity- 
noninvariant potentials gives values for p(0) which lie be- 
tween and 1. However, if one considers the sum of p(0) 
and ,5(0) in the uncoupled case, one obtains an even inte- 
ger. For a parity-invariant potential, i.e., V(—x) = V(x), 
p and p are equal and the zero-energy value of p will al- 
ways be an integer. 



III. SEGMENTED POTENTIALS AND 
FACTORIZATION OF THE S MATRIX 

For the Schrodinger equation without coupling it is 
well known that the reflection and transmission ampli- 
tudes satisfy a factorization formula. That is, if the po- 
tential is subdivided into a number of sections, then the 
total transmission and reflection amplitudes for the sys- 
tem can be expressed in terms of the amplitudes for each 
of the truncated pieces of the potential. The recent proof 
of Aktosun (22) may be generalized immediately to the 
case of N coupled equations. 

Following Aktosun, then, we subdivide the real 
line into J pieces. The boundaries of the segments 
are denoted by Xi, i=0, . . . , J, with —R — xq<xi< 
. . . <xj-i<x,j — R. The potential may then be writ- 
ten as a sum of truncated potentials as follows 



where 



.7-1 

V(x) = J2 V i(x), 



V(x) = { x j< x<x ]+i 
3 1 0, otherwise. 



(3.1) 



(3.2) 



The single indices on the truncated potentials should not 
be confused with the implicit double indices which label 
the various elements of the potential matrix. 
For a given j, then, let us define the matrix 
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-\k) -Tj\k)p 3 {k) 



Pj(kK\k) (rUk) 



(3.3) 



where Tj(k),Tj(k), pj(k),pj(k) represent the various am- 
plitude matrices for the truncated potentials with the 
usual boundary conditions. The amplitudes for the orig- 
inal potential are similarly arranged into a matrix, 



(*) -r'\k)p{k) \ 



m ~\p{k)r-\k) (ft(fc))- 1 J' 



(3.4) 



and the factorization formula is then simply given by 



j-i 



A(fc)= n^( fc )' 



(3.5) 



the form V(x) = v(x)M, where v(x) is a real- valued func- 
tion of x and M is a constant symmetric matrix. Such 
potentials have been used previously in various applica- 
tions (sec, for example, Ref. E6|). It is easy to prove 
that when the potential is diagonalizable by a constant 
orthogonal matrix, then all of the amplitude matrices (as 
well as their inverses and hcrmitian conjugates) commute 
with each other. 

If we use potentials of this type to construct finitely pe- 
riodic potentials with nonovcrlapping subpotentials, the 
analysis of Rozman et al. [ p3|p4] | follows in the same way 
for the matrix potential problem and the expressions for 
the amplitude matrices are straight forward generaliza- 
tions of their results. 



where the factors on the right side of the equation are 
ordered so that factors with lower subscripts occur to 
the l eft o f the ones with higher subscripts. The proof of 
Eq. (3.5) is completely analogous to that advanced by 
Aktosun in the N = 1 case and so we shall not review it 
here. The only added complication is the fact that the 
various amplitude matrices do not generally commute, 
but this has been properly accounted for in t he definitions 
of the A matrices. The utility of Eq. (3.5) will become 
apparent below when we use it to derive the amplitudes 
for scattering from two different delta-function matrix 
potentials in terms of the amplitudes for scattering from 
each of them separately. 

This approach effectively factorizes the S matrix in the 
sense that if the S matrix of the j th potential segment, 



Tj(k) 
Pj(k) 



Pj(k) 

Tj(k) 



(3.6) 



is known, then the scattering amplitudes of the potential 
segm ent are determined, and from them Ay(fe)- Using 
Eq. ( |3~5| ) we can obtain A(fc) for the whole potential, 
and this allows us to solve for the scattering amplitudes 
and the S matrix of the whole potential. 

Another generalization of the uncoupled to the cou- 
pled problem involves the finitely periodic potentials, re- 
in 



cently discussed by a number of authors [p|,[23| |2 
the derivation of the factorization formula, Eq. (3.5) 



have to be careful in the ordering of the products such as 
pr^ 1 . This non-commutativity of the amplitude matrices 
would typically prevent us from generalizing the closed- 
form solutions of the finitely periodic potentials. There 
are however classes of potentials for which the various 
amplitude matrices do commute with each other and for 
which the results for no coupling can be generalized. 

An example of such a class of potentials consists of 
those potentials which can be expressed as 



V{x) = Udia,g[vi(x),v 2 (x), . . .,v N (x)]U' 1 



(3.7) 



where U is a constant ( rea l) orthogonal matrix. Note that 
the potentials of Eq. (3.7) include as a subclass those of 



IV. POTENTIAL MODELS 

Below we consider some potential models which lend 
themselves to solutions in closed form. These models help 
to elucidate some of the results obtained in the previous 
sections. 



A. Constant potential matrix 

An example of a potential for which solutions can be 
obtained in closed form is the square-well or square- 
barrier potential matrix for which 



V(x) 



Vq for a < x < b where a > —R and b < R 
otherwise, 

(4.1) 



where Vq is a real s ymm etric N x N matrix. The 
Schrodinger equation (2.1) is equivalent to a first-order 
differential equation of the matrix function W{x) — 
W&,x) ofEq. ©.i.e., 



W'(x) = F(x)W(x) 



where 



F(x) 



V(x) 



1 





(4.2) 



(4.3) 



with the boundary condition W(—R) = I. As we saw ear- 
lier the function W(x) provides the advantage of giving 
solutions th at a rc linearly independent at k = 0. In ad- 
dition, Eq. (4.2) gives us an initial value problem, rather 
than the two-point boundary condition problem of the 
original Schrodinger equation. Solving for the scatter- 
ing amplitudes numerically is simpler for the initial value 
problem. An alternative to this approach is the variable 
amplitude formulation which also casts the problem into 
a system of first order differential equations with an ini- 
tial value condition |T^JI^]. In principle, Eq. (4.2) can be 
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used to solve the Schrodinger equation for any arbitrary 
potential matrix. 

For the con stant potential matrix the differential equa- 
tion Eq. (4.2) can be integrated starting at x = —R over 
the three regions (— R, a), (a, b) and (b, R) in turn p7| . 
The result is 



W{R) = 



cos k(R-b) k-T-awk^R-b) 
-ksmk(R-b) cos k(R — b) 

cosh If (6 — a) K" 1 sinh if (6 — a) 
K sinh K (b — a) cosh K(b— a) 

cosk(a + R) k^ 1 sink(a + R) 
-k sin k(a + R) cos k(a + R) 



(4.4) 



where K 2 — Vq — k 2 . To simplify matters, but still to al- 
low us to study a model with a potential function without 
definite parity, we set a = — R, so that we obtain explicit 
forms for the wave functions at x = R. 



<P(k,R) 
<P'(k,R) 
X (k,R) 
x'(k,R) 



cos k(R-b) cosh K{R + b) 

+k- x sink(R - b)K sinh K(R + b) 
-ksbxk(R- b) cosh K(R + b) 

+ cos k{R~ b)K sinh K{R + b) 
cos k(R - b)K- x sinh K(R + b) 

+k- 1 sin k(R - b) cosh K(R + b) 
-k sin k(R - b)K- x sinh K{R + b) 

+ cos k(R~b) cosh K (R + b). 



(4.5) 



These expressions can be inserted in the equations for 
to ( p. 51 ), to obtain the scattering 



(2, 



P,P,t,t, Eqs. 
amplitudes. 

Note that K 2 = Vq — k 2 is a real symmetric matrix and 
may therefore be diagonalized by an orthogonal transfor- 
mation U, giving K 2 D = UK 2 U T . The diagonal matrix 
Kb has the square root of the diagonal elements of K D 
along its diagonal. Thus K = U T K]jU, which is a sym- 
metric matrix. Consequently, <p,<f>', X>x' are symmetric 
matrices, which leads to r = r . Furthermore, since each 
of the wave-function matrices or their derivatives at R 
is a power series (or polynomial) of the matrix K, the 
wave-function matrices commute. It is not difficult to 
show that in the case of b = R, i.e., when the potential 
function has even parity, p(k) = p(k). 

The threshold behaviour of the transition amplitudes 
for a potential lacking specific parity can be studied ex- 
plicitly with this model. Since Levinson's theorem in- 
volves the trace of the amplitudes at zero energy, we 
need to consider the diagonalized forms of the amplitudes 
only. We diagonalize each of the wave-function matrices 
of Eq. (4.5) using the same orthogonal matrix U for each, 



and we denote the diagonal matrices at zero energy as 



<P' D (0,R) 
Xd(0,R) 



diag(pi,...,pAr) 
diag(Ai, ...,Xn) 
diag(xi,...,a;Ar) 
diag^i,...,^). 



(4.6) 



When there is no half-bound st ate, w e obt ain, b y insert- 
ing these expressions into Eqs. ( 2.48| ) and ( 2.51 ), 



p(0) = p D (0) = -1 and r(0)=TD(0) = O. (4.7) 

When there are n half-bound states and det 0'(O, R) = 0, 
we write <f>' D (0, R) = diag(0, . . . , 0, A n +i, . . . , Xn)- Using 
Eq. ( |2.13D we find that 



<£ f (k, x)x(k, x) - (j>\k, x)x'(k, x) 



(4.8) 



In general the matrices <p and x are rea l an d for the con- 
stant potential matrix they are symmetric as well. Thus 

</>'(0, J R) X (0,i?)-^(0,i 1 , )x , (0, J R) = 
#j(0, R)xd(0, R) - MO, R)x' D (0, R) = -1. (4.9) 

From this relation it follows that x[ = 1/pi for i = 
1, . . . , n, and furthermore 



p D (0) = diag 
Pd(0) = diag 

t d (P) = Td(P) 
= diag 



1 



pi 



i 



p 2 



i+Pl 



(4.10) 



Pi 



■Pi 



l+p 2 



2pi 



2Pr, 



1+pf 



1 



Pn 



,0,...,0 



-1 
(4.11) 

(4.12) 



Clear ly the relation (2.83) is satisfied by Eqs. (4.10) and 
(4.11). For the parity invariant potential obtained by 
setting b = R in Eq. ([D| , 0(0, R) = x'(0,i?). Hence 
Pi = x\ for i = 1, . . . , N, and it follows that p\ = 1. Such 
a potential therefore yields transition amplitudes of the 
form 

p D (0) = p D {0) =diag(0,...,0,-l,...,-l) and 
T D (0) = diag(±l,...,±l,0,...,0). (4.13) 



When the i diagonal element of tf>E>(0,x) is an even 
(odd) function, then the i th diagonal element of T£>(0) 
will have a plus (minus) sign with the one. The converse 
is not necessarily true. 

Consider the special case of N = 1. For the par- 
ity invariant potential with a half-bound state, one has 
p(0) = and t(0) = ±1. The plus sign corresponds to 
4>(0, x) being an even solution and the negative sign to 
it being an odd solution. When there is no half-bound 
state, then p(Q) = —1 and r(0) = 0. In the case of a 
potential without definite parity, p(0) = p(0) = — 1 and 
t(0) = r(0) = when the potential does not support 
a half-bound state. When there is a half-bound state, 
r(0) and r(0) are not equal to zero, nor are the p's equal 
to —1. However, the sum of the p's is an integer, i.e., 
p(0) + p(0) = 0. These results, which are clearly valid 
for the square-well potential, are actually valid for any 
N = 1 potential function. 







B. Models involving delta- function potentials 

We now consider two examples involving delta func- 
tions for which results can be obtained in closed form. 
The results exhibit qualitative features which are also 
found in much more complicated examples. First we ex- 
amine the case of a single delta-function matrix potential 
positioned at the origin. Then we will use the factoriza- 
tion formula derived earlier to look at the case for which 
there are two delta-function matrices symmetrically po- 
sitioned on both sides of the origin. 



1. Delta function at the origin 



We write the Schrodinger equation for this case as 



— +8(x)\) V = fc>, 



(4.14) 



where A is an N x N symmetric matrix, and tp can be 
taken to be either a column vector solution or a solu- 
tion matrix. The former approach will be used when we 
consider bound states and the latter when we examine 
scattering solutions. 

First consider the scattering solutions. Since the po- 
tential has even parity, we immediately have the result 
that p = p and r = r. Thus we need only consider the 
solution with the incident wave from the left, 



tp(k,x) 



le l 



„ikx 



pe 



~ lkx , x<0 
x>0. 



(4.15) 



Here we see the utility of working directly with a ma- 
trix of column eigenvectors (as opposed to working with 
individual column vectors); p and r may be solved for 
directly in terms of matrix operations. The scattering 
amplitudes are 



p(k) = (2ik-X)~ 1 X, 

r(fc) = 1 + p(k) = 2ik(2ik - A)" 1 . 



(4.16) 
(4.17) 



Of particular interest to us, due to its connection with the 
version of Levinson's theorem given in Eq. (2.43), is the 
quantity Tr[p(0) + p(0)). If A" 1 exists, then p(0) = -1 
and Tr[p(0) + p(0)] = 2Tr[p(0)] = -2N, as expected. If A 
is not invertible, however, we must be a bit more careful. 

In order to determine the significance of the nonin- 
vertibility of A, consider the bound-state solutions of Eq. 
( 4.14 ). Setting a 2 — —k 2 and insisting that a > 0, we 
find that the column eigenvector for the bound state is 



Ae c 
Ae 



x , x<0 
ax , x>0, 



(4.18) 



whe re A is a normalized column matrix. By integrating 



Eq. (4.14) over an infinitesimal interval including the ori- 
gin, we obtain an expression between the derivatives of 



^b{ot,x) on both sides of the origin, which leads to the 
relation 

{2a + X)A = 0. (4.19) 

In order to avoid the trivial solution, we demand that 

det(2a + A) = 0. (4.20) 

The non-negative values of a which solve the above equa- 
tion define the bound-state energies. Clearly there is at 
least one half-bound state if det A = 0. 

Returning to the scattering problem, we find that the 
easiest way to proceed is to first diagonalize the matrix A. 
Since A is real and symmetric, the diagonalization can be 
accomplished by using an orthogonal matrix U, so that 

X D = UXU~\ (4.21) 

where A^> is diagonal and U T — U~ l . If we now define 

ijj D = UipU- 1 , (4.22) 



we see that Eq. (4.14) may be rewritten as 
d 2 



dx 



2 + S(x)X D ) ip D = k 2 ip D . 



(4.23) 



The orthogonal transformation similarly transforms the 
boundary conditions, Eq. ( |4.15| ), to 



tp D {k,x) 



le lkx + p D (k)e~ ikx , 
T D (k)e lkx , 



x<0 
x > 0. 



(4.24) 



We see here another advantage of working directly with 
square matrices. If we had been working with column 
vector wave functions, the transformed wave functions 
would have been given by U^S, so that the normaliza- 
tion of the incoming wave would in general have been 
changed. Working with N x N wave-function matrices 
gives the above result that the form of the boundary con- 
ditions is unchanged under the transformation. In fact 
the transformed wave function is itself a diagonal matrix, 
and we essentially have N decoupled copies of the prob- 
lem with no coupling, with (possibly) different potential 
strengths.^] 

Suppose now that det A = 0. Then it follows that A 
has at least one zero eigenvalue. Let us again suppose 
that there are in fact n zero eigenvalues, so that 



X D = diag(0, 



, 0, A T . 



. Xn), 



(4.25) 



2 Note that this same trick can be employed any time the 
potential is of the form V(x) = v(x)M, where M is a real 
symmetric matrix. Diagonalizing M gives iV decoupled sys- 
tems with potentials Vi(x) — m,iv(x), i = 1, . . . , N , where the 
mi are the eigenvalues of the matrix M. 
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where the Aj, i — n + 1, ...,N, are the remaining 
(nonzero) eigenvalues. Then the diagonalized reflection 
and transmission amplitude matrices are given by 



p D {k) = diag 0, . . . ,0 



A n +i 



A 



JV 



2ik — A n _)_i ' 2ik — A^r 

(4.26) 



and 



T D (k) = diag 



so that 



2ik 



2ik 



2ik — A rl +i ' ' 2ik — Ajy / ' 

(4.27) 



Tr[p(0) + p(0)] = 2Tr[p D (0)] = -2(JV - n). (4.28) 

Thus we see in this example how the trace of p(Q) + p(0) 
keeps track of the number of half-bound states in the sys- 
tem. In fact it is easy to verify that Levinson's theorem 
holds for the coupled system, since it holds separately for 
each decoupled equation of the diagonalized problem]^] 
It is instructive to consider the relation 

2ikTp 1 (k) = diag(2ifc, . . . , 2ik, 2ik — A n +i, . . . , 2ik — \n), 

(4.29) 



which follows from Eq. (4.27). It demonstrates fo r thi s 
model that in the limit as k — > the expression Eq. ( 2.65 ) 
is real, as expected. 



2. Potential with two delta functions 

We now turn to a slightly more complicated example, 
in which there are two delta-function matrix potentials, 
one at x = a and the other at a; = —a. The N = 1 version 
of this model was studied by Senn 0] . The Schrodinger 
equation for this case is given by 



— -^+d(x + a)X + S(x-a)X)tp = k^ijj, (4.30) 
dx z J 

with boundary conditions 



for the wave incident from the left and 

Mk x) - {T(k)e- tkx , x<-a 
nK ' X) -\le- ikx +p(k)e ikx , x>a, 



(4.32) 



for the wave incident from the right. Rather than solve 
the Schrodinger equation again, we may now simply sub- 
stitute the results of the previous section into the factor- 
ization formula, Eq. ( |3.5|) f] An evaluation of the resulting 
expressions yields 

p(k) = (\e- Uka + (2ik + \)(2ik - A)~ ^ 

r _1 (fc, a; A, A)A _1 , (4.33) 
r(fc) =-Ak 2 e- 2lka {2ik-\)- 1 T- 1 {k : a\\,\)\- 1 : (4.34) 
p(k) = (2ik - A) -1 r _1 (A;, a; A, A) 

f(2ik - \)\- 1 \e- itka + 2ik + \\ , (4.35) 
f(fc)=r T (fc), (4.36) 
where 

T(fc, a; A, A) = (2ik - A)A- 1 e^ 2lfca - (2ik - A) _1 Ae 2ifeo . 

(4.37) 

Let us assume for the moment that both A -1 and A -1 
exist so that the above expressions are well defined. (For 
k > 0, it is actually sufficient that only one or the other 
exists - it is possible to rewrite the expressions so that 
they contain only A -1 and not A -1 .) Performing a Taylor 
expansion of p and p for small k, we find that in the 
typical case p(0) = p(0) = -1, so that Tr[p(0) + p(0)] = 
—2N, as expected. The atypical case is defined by the 
condition det(A _1 + A -1 + 2a) = 0, which, as we shall 
see, is also the condition for a half-bound state. 

Let us then work out the bound-state condition. This 
may be done in a manner similar to that for the single 
delta-function case to obtain 



det T(ia, a; A, A) = 0. 



(4.38) 



Solutions of Eq. ( 4.38 ) with a > correspond to bound 
states. As a — > 0, Eq. (4.38) yields the half-bound-state 
condition, 



det(A~ J + A" 1 +2a) = 0. 



(4.39) 



Alternatively, if we employ the wave functions <fi and x 
of Sec. II C for the model potential and use Eq. (2.58|) as 



the condition for the bound state, we obtain the equation 
det ( (A + 2a) (A + 2a) - AAe~ 4cM /2a) = 0, (4.40) 
which in the limit as a approaches zero reduces to 



3 The proof follows on noting that the determinant of the S 
matrix is unchanged under the transformation which diago- 
nalizes p and r. 



4 Some care must be taken with the reflection amplitude ma- 
trices, for they acquire a phase when the potential is trans- 
lated. The transmission amplitude matrices are, however, 
unchanged. 
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det(A + A + 2aAA) = 0. 



(4.41) 



This equation is preferred over Eq. ( 4.39| ) since it is not 
artificially singular when one of the inverse matrices does 
not exist. 

Let us now consider an explicit example with N = 
2. Since one of the two matrices A or A may always be 
diagonalized by an orthogonal transformation, we will 
let A be diagonal right from the start. An example which 
gives a half-bound state for a = 1 is one for which 



A 



-1 



A = 



(4.42) 



Fig. [l] shows a parametric plot of pu(k) as a func- 
tion of k in the complex plane for the cases a = 0.95, 
a = 1.00, and a = 1.05. In the two typical cases 
(a = 0.95, 1.05), pn(0) = —1, while for the atypical case 
(a = 1), pn(0) = 0.777. . .. This is then the analog of 
Senn's "threshold anomaly" for the generalized matrix 
version of his model [Q . Examination of the other diago- 
nal reflection amplitudes yields the expected result that 
Tr[p(0) + p(0)] is equal to —4 in the typical case and —2 
in the case where one half-bound state exists. 



0.5 



-0.5 




-0.5 





Re( p n ) 



0.5 



FIG. 1. Plot of the Im(pn(fc)) versus the Re(pn(fc)) for 
k — to k = 5 for double delta-function matrix potential. 
For the a — 1.05 case the curve reverses the direction of travel 
around the origin when k ~ 1.1. 

The behaviour of pn{k) as a function of k may strike 
the reader as being somewhat peculiar: for sufficiently 
large k as k increases, pu{k) traces out a never-ending 
counter-clockwise spiral towards the origin. A plot of 
the argument of pn(0) for the three cases would show 
that the phase shifts are not bounded - they keep on in- 
creasing to infinity. This peculiar feature does not exist 



when there is no coupling (the phase shifts are bounded 
due to Levinson's theorem), but is a rather generic fea- 
ture of N > 1 models. The important thing to bear in 
mind when N > 1 is that the phase shift which obeys 
Levinson's theorem is defined as being proportional to 
the logarithm of the determinant of the S matrix. This 
phase shift can in general be a nontrivial function of the 
"physical" phase shifts associated with the scattering am- 
plitudes. 




FIG. 2. The determinant function for the double 
delta-function potential graphed as a function of a. 

Fig. |^ shows a plot of 

f(a, a; A, A) = det ([(A + 2a)(A + 2a) - AAe" 4Qa ]/2a) , 

(4.43) 

as a function of a for a « 1. The inset in this figure 
shows an expanded view of the function near the origin 
for the three cases a = 0.95, 1.00, and 1.05. In addition to 
the two regular bound states that all three cases possess 
(near a — 0.5164 and a = 3.3508), the a = 1.05 case has 
an extra bound state near a = 0.0259, and the a = 1.00 
case has a new bound state just emerging at a = 0. 

Finally, Fig. 3 shows a plot of the "Levinson's theo- 
rem" phase shift as a function of k for the three cases. 
Clearly this phase shift is well-behaved and is bounded. 
As the potential "strength" is adjusted so that the sys- 
tem goes through a half-bound state the phase shift at 
k = jumps by tt in two increments of 7r/2. 
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0.5 



"1 1 
a = 0.95 

a=1.00 

a=1.05 



0.2 



0.4 



0.6 



0.8 



FIG. 3. The phase of the S matrix (divided by n) of the 
double delta-function potential as a function of k. For k larger 
than shown on the graph the three curves remain close to one 
another and approach zero as k ~ > oo. 



V. DISCUSSION 

In this section we make a few observations. The prob- 
lem of one-dimensional coupled-equation scattering us- 
ing a representation of wave functions which have incom- 
ing waves from the left or the right is readily solvable. 
Despite the advantages of a "partial wave" representa- 
tion QJ^] or a parity-eigenstate representation |^|| for 
parity-invariant potential functions, our analysis (which 
is valid for any potential) is quite manageable. 

The use of wave function matrices (see Refs. [ fl6|]28| 
for three-dimensional scattering and also |l3|,|l4| for onc- 
dimensional scattering) rather than column- vector wave 
functions, leads to simplified notation for a number of 
relations, e.g., the closure relation, Eq. ( |2.3C| ). One also 
finds that performing a unitary transformation on the 
scattering wave function matrix does not alter the nor- 
malization of the incoming waves, whereas it does for 
column wave functions. The introduction of the real ma- 
trix wave function solutions <f> and \ nas two distinct 
advantages. In the first place the Schrodinger equation 
for the scattering problem can be reduced to a system of 
first-order differenti al e quations with one-point bound- 
ary conditions, Eq. (4.2). The scattering amplitudes (and 
the S matrix) are algebraic expressions of these functions 
evaluated at R. Furthermore, unlike the solution matri- 
ces ip and "0, matrices (f> arid \ have linearly independent 
columns at threshold and consequently are convenient for 
investigating threshold behaviour. 



Our starting point with the wave functions ip and 0, 
which gives the definition of the reflection amplitudes p 
and p, yields a generalized and simplified understand- 
ing of threshold behaviour. Whereas previous work j|j7]] 
indicates that for parity-noninvariant potentials the re- 
flection amplitude at threshold can have noninteger val- 
ues, unlike that for parity-invariant potentials, we find 
that Tr|p( fc) + p(k)] at threshold is always an integer 
(see Eq. (f^83j)). The results of de Bianchi @, how- 
ever, already imply such a relation, as well as noninte- 
ger reflection amplitudes at threshold, for some parity- 
noninvariant potentials with no coupling. 

Finally, the phases of the reflection and transmission 
amplitudes of the coupled system are not simple func- 
tions of k, as is the case for uncoupled scattering, for 
which these phases satisfy the appropriate form of Levin- 
son's theorem. The phase r\ of the S matrix, which ap- 
pears in Levinson's theorem, is in general a nontrivial 
function of the phases of the scattering amplitudes. Fur- 
thermore, the function rj(k) is bounded, unlike the phases 
of the scattering amplitudes which are not necessarily 
bounded. 
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APPENDIX A: SOME PROPERTIES OF 
SOLUTIONS OF THE SCHRODINGER 
EQUATION 

In this appendix we determine the condition for linear 
independence of the solutions of the Schrodinger equa- 
tion. Consider the matrix Schrodinger equation 



0+F(^ = fcV, 



(Al) 



where ip is the N x N solution matrix whose columns 
are solutions to Eq. (2.1). Suppose we have two such 
solution matrices, / and /. We define a 2N x 2N matrix 
functional 



W(fJ) 



f f 

r f 



(A2) 



in which the prime indicates the der ivative with respect 
to x. Since / and / satisfy Eq. ( |Al|) the matrix W is a 
solution of the matrix equation, 



W = FW, 



(A3) 
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where 



F 





V-k 2 



(A4) 



Lemma 1 (det W)' — for all x G (—00,00). 

Proof: Let us write / and / in terms of their TV- 
component row vectors: 



( fA 

h 



and / 



V/W 

We then obtain [p7l page 96] 



( ~ h \ 

ft 

\f N ) 



(A5) 



(det W)' = det 



det 



ft 


ft 


h 


ft 


fN 


In 


ft 
ft 


ft 
ft 


fff 


/at 


ft 


ft 


ft 


ft 


fN 


In 


ft' 


ft' 


ft 


ft 



det 



det 



\ In In J 



zero because they have two equivalent rows. In order to 
show that the remaining terms are also zero, we write the 



Jl 


M 




ft 


ft 




fN 
ft 


fN 
ft 




ft 


ft 




\f'N 


In) 






M 




ft 


ft 




fN 


In 




ft 


ft 




ft 


ft 




\f N 


In) 






(A6) 


1 the 


right 


are 



Schrodinger equation, Eq. (Al), as 



ip"(x) = M(x)tp(x), 



(A7) 



where M(x) is an N x N matrix. Furthermore we write 
M[x) in terms of row vectors: 



M{x) = 



(A8) 



/ mi(x) \ 
m 2 (x) 

\ m N (x) J 

Then the matrix in the i second-derivative term of 



(ft h\ 

fN fN 
ft ft 



ft' ft' 
\f'w fN J 



1 







\ 




1 










1 




rrii{x) 




(J 








1 










' 1 / 



(W) 



(A9) 



The det erminant of the matrix on the right side of 
Eq. ( |A9| ) is zero, and consequently the (det W)' = 0. 
The det W is a constant function of x. 

Lemma 2 The solutions contained in the columns of f 
and f are linearly independent if and only if det W =/= 
for all x G (—00, 00). 

Proof: First suppose that det W ^ 0. We consider a lin- 
ear combination of solutions which is equal to the trivial 
solution, 



/h + /h = o for all x, 



(A10) 



where h and h are ./V-component column vectors of con- 
stants and o is the ./V-component zero column vector. A 
similar relation holds for the derivatives of / and /, so 
that 



where 



Wc = o, 



(All) 



(A12) 



Eq. (A6) may be written as 



and c and o are now 2iV-component column vectors. If 
det W 7^ for some x, which according to the previous 
result means it is nonzero for all x, c = o and the column 
solutions contained in / and / are linearly independent. 

Suppose now that det W = 0. If det W — for some 
x = xq, then the system of linear equations W(xq)c = o 
has a nontrivial solution c. We form a col umn solution 
of the system of differential equations, Eq. ( |A3| ), w(x) — 
W(x)c wh ich vanishes at xq- This is the trivial solution 
of Eq. (A3); w(x) = o for all x. It follows that the column 
solutions contained in / and / are linearly dependent. 



APPENDIX B: ANALYTIC PROPERTIES OF 
THE SOLUTION MATRICES AND X 

We co nsid er solution matrices <j){k,x) and x(fc,x) 
of Eq. (Al) with boundary conditions (f>(k,—R) = 
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X'(k,-R) = 1 and (f>'(k,-R) = x(k,-R) = 0, where 
R is the range of the potential. According to a theorem 
of Poincare an ordinary differential equation containing 
an entire function of some parameter has solutions which 
are entire functions of the parameter provided these so- 
lutions have boundary conditions which are independent 
of the parameter. We will show that </> and x are entire 
functions of k, following a similar derivation for partial- 
wave solutions in three-dimensional scattering p8| , p9| . 

It is straightforward to verify that the matrix functions 
4> and x with the given boundary conditions are solutions 
of integral equations of the Volterra type, 



<t>{k,x) 



X{k,x) = 



1 cosfc(a; + R) 

. sin k(x — x') 
ax 

sin k(x + R) 
k 

f sin k{x — x') 



dx 



V(x')(f>(k,x') (Bl) 



V(x') X (k,x'). (B2) 



-R 



In order to show that each element of the solution matri- 



ces is an entire function of fc, we rewrite Eq. (Bl) in the 
form 

4>(k, x) = 1 cos k(x + R) 



dx' dt cos/rt V(x')4>(k,x') (B3) 

R JO 



We solve Eq. (B3) by successive approximations of the 
form 



4>(k,x) =J2^ s Hk,x), 



(B4) 



s=0 



where 

/X r-X — x' 

dx' \ dt cosfci VV)^" 1 ^^'), 
-R JO 

for s > 1 and </> (0) (k, x) = 1 cos k(x + R). (B5) 
Thus 



\^(k,x)\< 



dx' 

R JO 



dt | cos kt\ 



and 



\4?(k,x)\ =Sij\ cos k(x + R)\. 



s > 1 (B6) 



(B7) 



Denoting Sfc for the imaginary part of k and using the 
relation | cos fct| < cosh3fc£ for t real, we obtain upon 
iteration 



i<4 s) (m)I<Ei41'(m)i 

i'j> 

sinh(23/fci?) 



< 



/x rti 
dxi / dx2 ■ ■ ■ 
-R J-R. 



x s -i 
R 



dx s E \ v i'lM\\ v hhfa)\''-\Vi.j>(x,)\. 

i',j',h,—,la 

(B8) 



Since the integrand is a symmetric function under the 
interchange of any pair (xi, Xj), 



\^(k,x)\< 



E 

i',j',h,...,l, 

Let 



sinh(2Qfci?) \ ' 
dxi {Vi'i^xi)] ■ 



cosh(23A; J R)^ T 



-R 



dx s \V lmj ,(x,)\. (B9) 



M = max / dx' \V tj (x')\ < oo 
l 'i J-r 

Then 

l4 s) (M)l< 



(BIO) 



sinh(23fci?) \ s , . _ , nN N s+2 , rs 
1 cosh(23fci?) — M S 



3fc 



Thus the series J2 S "Ay converges absolutely and uni- 
formly for x £ [-R, R] and for every region in the com- 

plex k plane. To determine the existence of — — (k,x), we 

ok 



(Bll) 



differentiate Eq. (BS) with respect to k, 

d<t>, 



dk 



(k,x) = -{x + R)lsink(x + R) 



dx' / dt t sin kt V(x')(j)(k, x') 
R Jo 



I*X px — x pti 

+ dx' dt coskt V(x')-£(k,x'). (B12) 

Jr Jo ok 



When Eq. (B4) is differentiated with respect to k it yields 

(B13) 



in- 



to \ V^ W r* i 



s=0 



where now we have 



dk 

with 
dk 



-(k,x) 



dx' 



R 



dt coskt V(x')- 



(8-1) 



dk 



<k,x') 
(B14) 



(k,x) = ~(x + R)lsmk(x + R) 



dx' dttshxkt V(x')(f>{k,x'). (B15) 

Jo 
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-{k,x) 



is bounded, 



and the convergence of series ( P13| ) follows in the same 
manner as that of (j)(k,x). Since 4>(k,x) and its deriva- 
tive with respect to k exist for all k, 4>{k,x) is an entire 
function of k. Similarly x ) can be shown to be an 
entire function of k. 

For real k the behaviour of 4>{k, x) and x(^) x ) as 
k becomes very large can be determined by iterating 
Eqs. @ and (jB2l) . Thus 



4>(k,x) ,~ 1 cosk(x + R) + 
1 f x 

- / dx'smk(x- x')V(x')co8k(x' + R) + 0(l/k 2 ) 
k J-r 



(B16) 



and 



X(k,x) 



-R 



sin k(x + R) 



dx' sin k(x - x')V{x') sin k{x' + R) + 0{l/k 3 ). 
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